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ABSTRACT 

Clustering of cosmic-ray sources affects the flux observed beyond the cutoff imposed 
by the cosmic microwave background and may be important in interpreting the 
AGASA, Fly's Eye, and HiRes data. The standard deviation, a, in the predicted 
number, N, of events above 10 20 eV is a/N = 0.9(ro/10 Mpc) 0,9 , where tq is the 
unknown scale length of the correlation function (r*o — 10 Mpc for field galaxies, 
Hq = 50 km s _1 Mpc -1 ). Future experiments will allow the determination of tq 
through the detection of anisotropics in arrival directions of ~ 10 20 eV cosmic-rays 
over angular scales of ~ tq/30 Mpc. 



1. Introduction 

The conventional astronomical picture for the origin of ultra-high energy cosmic rays, namely 
proton acceleration to high energy in extra-galactic objects, predicts a sharp suppression of the 
cosmic-ray flux beyond the Greissen-Zatsepin-Kuzmin (GZK) cutoff at ~ 5 x 10 19 eV (Greisen 
1966; Zatsepin & Kuzmin 1996), due to interaction of protons with photons of the cosmic 
microwave background. The absence of a GZK cutoff might suggest the presence of a new source 
of ultra-high energy cosmic rays, possibly related to the decay of exotic particles (see Cronin 1996; 
Hillas 1998 for reviews). 

The suppression of flux beyond the GZK cutoff is most often discussed assuming a uniform 
source distribution. However, the distribution of other astrophysical systems (e.g. galaxies, clusters 
of galaxies) is inhomogeneous on scales of tens of Mpc, comparable to the propagation distance 
of protons of energy > 5 x 10 19 eV. Therefore, significant deviations from the predictions for a 
uniform distribution may be expected beyond the GZK cutoff (Giller, Wdowczyk & Wolfendale 
1980; Hill & Schramm 1985; Waxman 1995). 

There are at least two possible approaches to evaluating the effects of source inhomogeneity 
on the energy spectrum and spatial direction of high energy cosmic rays. In the first approach, 
one assumes that the source density of ultra-high energy cosmic rays is proportional to (possibly 
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with some bias factor) the galaxy density in some particular survey of the distribution of relatively 
nearby galaxies. This approach has been used by Waxman, Fisher & Piran (1997); Giller, 
Wdowczyk & Wolfendale (1980);and Hill & Schramm (1985) and illustrates some of the principal 
features of source clustering. In the present paper, we adopt a different and complimentary 
approach. We use an analytic model that summarizes the clustering properties of the unknown 
source of ultra-high-energy cosmic rays by a single parameter, vq, the correlation length. 

The analytic model that we adopt has the advantage of generality. No specific source 
population has to be assumed. Also, the results of source clustering on the energy spectrum 
and on the spatial distribution of the highest energy cosmic rays can be summarized in terms of 
the single unknown correlation length of the source population, vq. Thus the correlation length 
provides a concise and simple characterization of measurable clustering effects on observational 
parameters for high energy cosmic rays. Future experiments that measure the arrival directions 
of a large number of high energy cosmic rays will determine r$ and permit a clear assessment of 
whether or not the departures from expectations based upon a homogeneous source distribution 
can plausibly be explained by any non-uniform population. 

The source correlation function, £(df), is defined by an average over the observable universe 
volume, (n(r)n(f+ dr)) = (n) 2 (l + £(dr)) where n is the source density and brackets denote 
volume average. Thus, the variance that we derive in the expected cosmic-ray number due to 
source clustering is, strictly speaking, the variance of the distribution of cosmic-ray number 
observed by observers randomly distributed over the universe. This is the best that one can do, 
but it is not exactly what we we want. We really want the variance in conceivable realizations 
of the universe in our vicinity. However, we show in §2 that the distance from our Galaxy over 
which sources contribute to the observed flux above 10 20 eV is ~ 40 Mpc (see Fig. 1). Since the 
clustering properties of astronomical objects (e.g. galaxies) within this volume are similar to that 
obtained by averaging over larger volumes (e.g. Peebles 1999), we expect our result to provide a 
reasonable approximation for the variance in the number of cosmic-rays observed at Earth due 
to different realizations of cosmic-ray source distributions that satisfy the universal clustering 
properties of the sources. 

One may argue that a more accurate estimate of the variance can be derived by the using 
known galaxy catalogs to model the effects of source inhomogeneity. However, in order to use 
such an approach one must choose a specific model to derive the probability of a given cosmic-ray 
source distribution under the constraints provided by a given galaxy catalog. The variance 
calculated in this approach would depend primarily on the resultant cosmic-ray source correlation 
function. While this correlation function would indeed depend on the clustering properties of 
the particular galaxies chosen, it would be strongly dependent on the assumed model. It should 
be emphasized here that one does not known which particular rare sub-population of the total 
observed galaxy population actually produces the highest energy cosmic rays. The local density of 
field galaxies is of order 0.1 Mpc -3 , whereas the local density of sources of UHE cosmic rays may 
be as low as 0.00002 Mpc -3 . (The lower limit to the source density is set by the requirement that 
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at least a few sources exist out to a distance of ~ 40 Mpc in order to account for the observed 
events above 10 20 eV.) Thus, results derived in this approach would not yield a more accurate 
estimate of the variance, but rather reflect the assumptions under which the cosmic-ray source 
distribution is constrained given some chosen galaxy catalog. Moreover, even after one makes 
a specific assumption about the functional relation between the density of sources of ultra-high 
energy cosmic rays and an observed sample of galaxies, one has to perform and analyze many 
Monte Carlo simulations with galaxy catalogs in order to determine the statistical effects on the 
observed characteristics of ultra- high energy cosmic rays. 

We believe therefore that the simplicity and conciseness of the analytic model approach 
justifies its application, as one of the possible approaches, to the evaluation of the effects of source 
inhomogeneity. 

In this paper we consider the implications of clustering of cosmic-rays sources, adopting the 
conventional picture of protons as the ultra-high energy(UHE) cosmic-rays. We assume that the 
correlation function between the sources of ultra-high energy cosmic rays has the same functional 
form as for galaxies, clusters of galaxies, and quasars, but we do not require — as was done in 
previous work (Waxman, Fisher & Piran 1997; Giller, Wdowczyk & Wolfendale 1980; Hill & 
Schramm 1985) — that the distribution of cosmic ray sources be the same as for nearby galaxies. 
As described above, the conditions that must be satisfied in order to produce ultra-high energy 
cosmic rays are so exceptional that the spatial distribution of sources of UHE cosmic rays could be 
very different than the distribution of average, nearby galaxies like the IRAS sources discussed in 
Waxman, Fisher & Piran (1997). Also, we are primarily concerned in this paper with the energy 
spectrum of UHE cosmic rays, whereas the IRAS galaxy distribution was used in Waxman, Fisher 
& Piran (1997) to discuss their possible angular distribution on the sky. 

Our principal result is that the standard deviation due to clustering, o" c i uster mg(-E'c), in the 
number of cosmic rays detected above energy E c is proportional to the number, iV smoo th(-Ec), 
predicted for a uniform source density. This behavior contrasts with the standard deviation due 

1/2 

to shot noise, which becomes proportional to N ' ^ ^he cons tant of proportionality between 
"clustering and -/V" smoot h depends upon energy and is calculated in Section 2. For E c ~ 10 20 eV, the 
constant of proportionality between <T c i us tering and iV smoot h is of order unity for plausible values 
of the correlation length, vq, of the correlation function of the sources of the UHE cosmic rays. 
Anisotropics in the source distribution on a scale of ~ (ro/100Mpc) should be detectable in 
the angular distribution of ~ 10 19 ' 7 eV cosmic rays with the high rates that will be observed 
in the HiRes (Corbato et al. 1992), the Auger (Cronin 1992; Watson 1993) and the Telescope 
Array (Teshima 1992) experiments. Unlike the predictions of many particle-physics explanations 
(for recent reviews see Berezinsky 1998; Bhattacharjee 1998) of the Fly's Eye (Bird et al. 1993, 
1994) and AGASA (Hayashida et al. 1994; Yoshida et al. 1995; Takeda et al. 1998) results, no 
characteristic dependence on Galactic coordinates is expected on the basis of this conventional 
extragalactic scenario. 
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The computed large value of o" c i uster ing reflects the fact that the universe, and hence the 
UHE source population, is inhomogeneous over distances comparable to the mean free path 
of UHE protons that move through the cosmic microwave background radiation. The ratio 
^clustering /-^smooth cannot be reduced by observing longer or doing different experiments. We only 
have one "nearby universe" and we do not know a priori the local clustering properties of the 
sources of UHE cosmic rays. 

The angular scale of anisotropies due to magnetic scattering of UHE cosmic rays from 
individual sources is smaller than what is expected from the clustering of the sources themselves 
unless the inter-galactic field is close to the maximum value allowed with the available data 
on Faraday rotation. The upper limit on the contribution of an inter-galactic magentic field 
to the Farady rotation of distant sources, RM< lrad/m 2 for sources at z = 2.5 (Vallee 
1990), implies an upper limit B < 10 -10 G on a uniform inter-galactic field, and an upper 
limit BX 1 / 2 < l(T 8 GMpc 1/2 on a field with correlation length A.[] A simple random walk 
calculation shows that the upper limit on BX 1 / 2 sets an upper limit to the magnetic deflections, 
9 S < 0.04((i/lMpc) 1 / 2 (£;/10 20 eV)- 1 for propagation distance <i.| Galactic magnetic fields 
contribute only relatively small deflections. 

We present our principal results in Section & We show in Section ^ that the large value 
of o" c i us tering /-^smooth is important for the interpretation of the existing AGASA and Fly's Eye 
experiments. We also list in Section || the effects of clustering that may be detectable with the 
HiRes, Auger project, and the Telescope Array experiment. 



The value we quote is a few times larger than quoted by Vallee, since the upper limit is inversly propotional 
to the electron density n e and we choose n e — 3 x lCT 7 cm -3 , corresponding to Qth 2 = 0.03, while Vallee assumed 
n e = 10 -6 cm~ 3 . Our upper limit is much smaller than that recently claimed by Farrar & Piran (2000), mainly due 
to the fact that they use older radio data (Kronberg & Simard-Normandin 1976) for which the RM upper limit for 
z = 2.5 sources is weaker, RM< 5rad/m 2 , and assume that the free electron density is only 0.3 of the baryon density. 

2 Following the suggestion of Kulsrud et al. (1997), that magnetic fields could be amplified by turbulence associated 
with the formation of large scale filaments and sheets, several authors have recently considered propagation of cosmic- 
rays in large scale inter-galactic magnetic field (e.g. Sigl et al. 1999, Farrar & Piran 2000). The upper limit on 
magnetic field coherent on 1 Mpc scale in such structures is smaller by a factor f 1 ^ 2 than the upper limit of 10 -8 G 
on a volume filling field, where / is the fraction of the volume occupied by large scale filaments and sheets. This is 
due to the fact that the upper limit is inversly proportional to the electron density, n e oc / , and propotional to 
the square root of the path length I oc f of light through magnetized plasma. While this reduction of the upper limit 
may be partly compensated by assuming strong negative evolution of the magnetic field with redhsift (e.g. Farrar & 
Piran 2000), the assumption that the inter-galactic magnetic field is confined to large-scale structires does not lead 
to increase in the upper limit on the deflection of UHE cosmic rays. 
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2. Results 



The number, N, of high-energy cosmic rays predicted on average above a threshold energy E c 

is 

/■oo r 

N smooth (E c ) = (N(E>E C )) = / dE d 3 rP(r,E;E c ) 

((dS/dE) r ^n(r))/47Tr 2 , (1) 

where P (f, E; E c ) is the probability that a proton created at r with energy E arrives at earth 
with an energy above threshold. The energy spectrum generated by the sources is dS/dE and 
the luminosity- weighted source density is n{r). We present calculations for dS/dE oc E~ 2 and 
dS/dE oc E~ 3 , which spans the range usually considered in the literature. The AGASA and Fly's 
Eye results above 10 19 eV can be shown to require a spectrum less steep than E~ 2 ' s (Waxman 
1995). 

It is convenient to rewrite Eq. (|l]) in terms of a function F (r), where 



(N(E>E C )) = / drF(r)(n), (2) 
Jo 

and F is the survival probability averaged over energies, 

roc 

F(r) = / dEP (r, E; E c ) (dS/dE) . (3) 
JE C 

The constant of proportionality between dS/dE and E~ rn can be chosen arbitrarily since the 
constant cancels out of the ratio of standard deviation to expected number, which is the primary 
quantity we calculate (see Eq. [Tl] ). For convenience, we fix the constant so that F(0) = 1.0; we 
also suppress the dependence of F on E c . 

For ultra-high energy cosmic rays (E > 5 x 10 19 eV), the probability P(r, E; E c ) is determined 
by the energy loss of protons due to pair and pion production in interaction with the microwave 
background. The energy loss in a single pair production interaction is small, of order m e /m p , 
and for protons at energy > 10 19 eV the characteristic energy loss time due to this interaction is 
comparable to the Hubble time, ~ 5 x 10 9 yr (Blumenthal 1970). Hence, pair production has only 
a small effect on P(r, E; E c ) for E c > 10 20 eV, where the proton life time due to pion production is 
< 3 x 10 8 yr. We calculated therefore P(r, E; E c ) taking into account pion production interactions 
only. We used the compilation of cross sections given in Hikasa et al. (1992), and assumed that 
the total cross section corresponds to the isotropic production of a single pion (this assumption 
has been shown to be adequate in Yoshida & Teshima 1993; Aharonian & Cronin 1994). The 
resulting probability distribution may be approximated well over the energy range of interest by a 
function of the form 

P (r, E; E c ) = exp [-a (E c ) r 2 exp (6 (E c ) /E)] . (4) 

We give below results obtained using a precise numerical evaluation, of P(r, E; E c ). The 
approximation (|j) provides a simple description of the functional dependence of P on r and E and 
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can be used above 8 x 10 19 eV to evaluate to an accuracy of 10% the average quantities considered 
in this paper. The appropriate values of a and b for E c /(10 20 eV) = 1,3, and 6 are, respectively, 
a/(l(T 4 Mpc- 2 ) = 1.4, 9.2, and 11, 6/(10 20 eV) = 2.4, 12, and 28. 

The function F(r) is shown in Fig. 1 for a spectrum dS/dE oc E~ 2 . The average distance 
from which a proton originates if it is observed at earth to have an energy in excess of E c , 

_ JdrrFjr) 

{r)= JdrF(r) ' (5j 

can be calculated analytically using Eq. (|) and Eq. (|). For an energy spectrum with 
dS/dE oc E~ 2 , 

VA^l [1 - exp(-6/2£ c )] 

A similar expression can be obtained for dS/dE oc E~ 3 . Inserting the appropriate values of a and 
b in Eq. (g), (r) = 31.2 Mpc for E c = 1 x 10 20 eV and (r) = 10.6 Mpc for E c = 3 x 10 20 eV (the 
values obtained using the numerical calculation of F(r) are 31.8 Mpc and 10.9 Mpc, respectively). 
The observed universe is inhomogeneous on these distance scales. 

The variance in the number of cosmic rays observed above an energy E c can be computed 
from the expression 

o 1 (E > E c ) = (N 2 (E > E c )) -(N(E> E c )) 2 . (7) 

Since galaxies, clusters of galaxies, and quasars are all clustered with correlation functions that 
have the same shape, it is natural to suppose that the sources of UHE cosmic rays are also 
clustered with a similar correlation function. Explicitly, we assume that 

(n (f ) n (f")> = (n) 2 (1 + £ (f - r") + S(r' - r")/(n}) , (8) 

where the correlation function £ is 



- ) 0) 

and the 5 function contribution represents shot-noise. The correlation length, ro, is not known 
for the sources of UHE cosmic rays, but it is about 10 Mpc for galaxies (Groth & Peebles 1977. 
Shectman et al. 1996, Tucker, Oemler, Kirshner et al. 1997) and about 40 Mpc for rich clusters 
of galaxies (Bahcall 1985; Bahcall & Soneira 1983; Peacock and West 1992) (for H = 50 km s" 1 
Mpc" 1 ). 

When the variance is calculated using Eq. (|7j), the only term that survives, in addition to shot 
noise, is proportional to the correlation function, i.e., 



a 2 (E > E c ) = J J dE'dE" (dS/dE') (8S/dE") 
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x („(r')){n (f")>J(f'-r") + (10) 

For simplicity, we will not display o" s hot-noise explicitly in what follows. After some algebra, we 
find 



(5/2) 1 / 2 r °- 9 



•^clustering 
-^smooth 

■jjoo ; oo ^ F[x)F{y) | (x + y) 0.2 _ {{x _ j]°- 5 



Jo° drF{r) 



(11) 



The ratio clustering /-^smooth can be evaluated numerically using Eq. (|3|) and Eq. (|ll|). For 
cosmic rays above 1.0 x 10 20 eV, we find 

^clustering Q 9 f r ° ~\ ° 9 /-^2) 



^smooth VlOMpC 

and for cosmic rays above 3.0 x 10 20 eV, 



^clustering 2 Q ( ^0 \ (13) 



^smooth V 10 M P C 

The results given in Eq. ( |T^ ) and Eq. (12) were calculated assuming a spectrum oc E~ 2 . If the 
spectrum is instead assumed to be oc E~ 3 , then the coefficient 0.9 in Eq. dl2| ) is increased to 1.0 
and the coefficient in Eq. (|13|) becomes 3.6. For cosmic rays above 3.0 x 10 20 eV, the dispersion is 
so large (see Eq. |l3| ) that it will be difficult to interpret statistically the observed number. 

Due to the rapid decrease with r of F(r) at large r, shown in Fig. (1), the integral ( |i"o|) 
is dominated by the contribution from small separations \f* — r"\. Neglecting, for example, the 
contribution from large separations |f" — r"\ > 20 Mpc, where the galaxy-galaxy correlation 
function drops faster than r reduces clustering /iVsmooth by 20% for E c = 10 20 eV and by 3% 
for E c = 3 x 10 20 eV. Thus our results for cr cmsiierm g /iV srnoot h are n °t sensitive to the r dependence 
of £(r) at large r. On the other hand, 0" c i us tering /iVsmooth is more sensitive to the r dependence at 
small r. For £(r) oc r~ 1,8 and ro = 10 Mpc, approximately half the contribution to the integral 
( |To| ) comes from separations |r" — f"\ < 3 Mpc for E c = 10 20 eV and |r" — r"\ < 1 Mpc for 
E c = 3 x 10 20 eV. The galaxy-galaxy correlation function, for example, is known to follow the 
form (|9|) down to r ~ 0.1 Mpc, while the cluster-cluster correlation function is not well measured 
for r < 10 Mpc. 

In principle, the contribution to the variance from shot-noise, c 2 hot _ noise , could be large due to 
the possibility of having a nearby source. However, an extremely close source would dominate the 
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all sky flux, which is not the case for the observed cosmic-ray events. If the average source density 
is much larger than the minimum value required to explain the observations of UHE cosmic rays, 
i.e. < n > >> 10~ 4 Mpc~ 3 , then the shot noise contribution is relatively small. The shot noise 
could be significant, of order AT smooth , if the local source density is of order 10~ 4 Mpc~ 3 . 



Are new particles required to explain the observed number of ultra-high energy cosmic rays? 
The number of events observed beyond 10 20 eV by AG AS A (Takeda et al. 1998), by the Fly's 
Eye (Bird et al. 1993, 1994), and by Yakutsk (Efimov et al. 1991) is, respectively, 6, 1, and 1. 
The average observed number of events, 2.7, is to be compared with conventional estimates of 
the number of events from a smooth model of between Ng^^ = 0.7 (for dS/dE oc E~ 5 ) (cf. 
Takeda et al. 1998) and A smooth = 2.2(for dS/dE oc E~ 2 ) (cf. Waxman 1995). The total standard 
deviation in the expected number of events beyond 10 20 eV is 



where the first term in the brackets is the calculated uncertainty due to clustering and the 
second term is due to Poisson statistics. For N smoo th = 2.2 and ro = 10 Mpc, ertotai = 2.4. Even 
for iVsmooth = 0.7, o" tota i = 1.0. No matter how one combines the experiments and the smooth 
models, there is not a 3<7 discrepancy with the conventional picture of protons as the source 
of ultra-high energy cosmic rays. If one considers, as many recent authors have done, only the 
AGASA data, then A^bserved = 6 and for dS/dE oc E~ 2 , iV smoot h = 2.2, and <7totai = ^-4, which is 
a 1.6a discrepancy. On the basis of the available evidence, we conclude that new particles are not 
required to explain the observed cosmic ray energy spectrum. 

However, a number of authors have interpreted the same data as suggesting the possibility 
that new particles are producing the highest energy cosmic rays (for a recent discussion and review 
of this point of view see Ellis 1999 and the references contained therein). One of the principal 
reasons for the difference in conclusions is the uncertainty over whether the extrapolation of the 
'conventional high energy component' should be made using dS/dE oc E~ 2 , as we have done, 
or whether the fall-off of the spectrum at energies just below the GZK cutoff is stronger, e. g., 
oc E~ 3 . This is just one of the important questions that will have to be answered in the future by 
a much larger data set. 

We also note that preliminary results from the HiRes experiment were recently presented in 
a talk at TAUP99 (Matthews 1999), reporting 7 events beyond 10 20 eV for an exposure similar 
to that of the Fly's Eye. It is difficult to decide how this result should be interpreted, since the 
discrepancy between HiRes and Fly's Eye results is present not only above 10 20 eV, but also at 
lower energy, where Fly's Eye, AGASA and Yakutsk experiments are in agreement: 13 events 
above 6 x 10 19 eV are reported in the preliminary HiRes analysis, while only 5 events at that 
energy range are reported by Fly's Eye. We therefore believe that unambiguous conclusions 
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Discussion and Predictions 



tftotai = 0.8 (ro/10 Mpc) L8 iV 2 mooth + N smooth 



(14) 
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based on the recent HiRes data can be drawn only after a complete analysis of the HiRes data is 
published (which would include, e.g., corrections due to realistic atmospheric conditions). 

One possible interpretation of the suggested excess of events beyond the GZK cutoff is that 
this may be a hint that ultra-high energy cosmic rays are produced by sources whose density in 
the nearby universe is somewhat higher than their average cosmological density. Of course, other 
interpretations are possible, including the hypothesis-amply represented in the literature-that 
new physics is involved and that the GZK cutoff is bypassed by this new physics. If the GZK 
cutoff is ultimately convincingly measured, then the most likely explanation for the small number 
of already observed UHE cosmic rays will be the effects of source clustering. 

Future measurements of high statistical significance with AGASA, HiRes, Auger and the 
Telescope Array are required to determine if source clustering affects significantly the observed 
number and angular distribution of UHE cosmic rays. The predicted effects of clustering depend 
upon the unknown correlation length r^. 

We summarize below the principal observational implications of assuming a correlation length 
between 10 Mpc (observed for galaxies) and 40 Mpc (observed for the rich clusters of galaxies). 

1) Clustering gives rise to a standard deviation that is proportional to iV smoot h (see Eq. ( |i~2| ) 
and Eq. ( |13| ) ). When more events beyond the GZK cutoff are available, clustering will be more 
important than shot noise in determining the shape of the energy spectrum of the highest energy 
cosmic rays. Even for the existing Fly's Eye and AGASA samples, the variance due to clustering 
may dominate shot noise (see Eq. ([Li])). 

2) The fractional calculated variance, (^clustering /-^smooth? increases with the cutoff energy, 
E c , which simply reflects the fact that the mean distance from which protons originate 
decreases with the observed energy (see Eq |6|) . The fractional variance is large at all energies 
beyond the cutoff. With the approximations of the present paper, (r) ~ 120 Mpc and 
^clustering/ ^smooth = 0.25 (r /10 Mpc) ' 9 for E c = 6x 10 19 eV [the exact value of er/iV sm0 oth ((?"}) 
is somewhat larger since at this energy pair production energy loss, neglected in our calculations, 
becomes important]. The fractional variance does not increase significantly as E c is increased 
above 3 x 10 20 eV, since the average distance at which a proton originates is approximately 
constant at higher energy (e.g., (r) = 9.5 Mpc for £ c = 6 x 10 20 eV, comparable to (r) = 10.9 Mpc 
for E c = 3 x 10 20 eV). 

3) Anisotropics should be observed on large angular scales ©o ~ foAorigim where r or i g i n decreases 
from ~ 200 Mpc to ~ 30 Mpc as E c increases from 0.4 x 10 20 eV(just below the GZK cutoff) to 
1 x 10 20 eV. A comparison of anisotropies observed below and above the GZK cutoff will be a 
crucial test of whether clustering is important. For ro ~ 10 Mpc, one can show that one will have 
to observe 10 2 (10 3 events) above 10 20 (4 x 10 19 ) eV to detect a 3a effect. Thus, the detection of 
anisotropies would require the large exposure of the Auger project. 

4) No characteristic dependence on Galactic coordinates is expected. This is in contrast to many 
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exotic particle-physics scenarios for the production of ultra-high energy cosmic rays in which 
a strong dependence on Galactic coordinates is predicted, in particular, a peaking toward the 
Galactic center. The Fly's Eye experiment reports an angular distribution peaked towards the 
Galactic disk for lower energy cosmic rays, which we interpret to be due to a different, more local 
set of sources than the sources of ultra-high energy cosmic rays. 

5) The suggested correlation of the directions of ultra-high energy cosmic rays with the directions 
of compact radio-loud quasars (Farrar &; Bierman 1998) should disappear as more events become 
available. The candidates suggested are all at distances > 10 3 Mpc, much beyond the allowed 
range of ~ 40 Mpc for ultra-high energy protons. 

We are grateful to an anonymous referee for constructive comments on the presentation in the 
initial manuscript. JNB is grateful to D. Eisenstein for instructive conversations about correlations 
functions and to NSF grant $J?HY95-13835 for financial support. Much of this work was done 
during a visit by JNB to the Weizmann Institute of Science. 
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r (Mpc) 

Fig. 1. — The survival probability averaged over energies. The plotted function, F(r), represents 
the probability that a proton created at a distance r will reach the earth with an energy in excess 
of the cutoff energy E c . Results are shown for E c = 1 x 10 20 eV and E c = 3 x 10 20 eV with an 
energy spectrum dS/dE oc E~ 2 . 



